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ABSTRACT 

We present measurements of 5-25 /.tm emission features of brightest cluster galaxies (BCGs) with 
strong optical emission lines in a sample of 9 cool-core clusters of galaxies observed with the Infrared 
Spectrograph on board the Spitzer Space Telescope. These systems provide a view of dusty molecular 
gas and star formation, surrounded by dense. X-ray emitting intracluster gas. Past work has shown 
that BCGs in cool-core clusters may host powerful radio sources, luminous optical emission line 
systems, and excess UV, while BCGs in other clusters never show this activity. In this sample, we 
detect polycyclic aromatic hydrocarbons (PAHs), extremely luminous, rotationally-excited molecular 
hydrogen line emission, forbidden line emission from ionized gas ([Ne II] and [Ne HI]), and infrared 
continuum emission from warm dust and cool stars. We show here that these BCGs exhibit more 
luminous forbidden neon and H2 rotational line emission than star-forming galaxies with similar total 
infrared luminosities, as well as somewhat higher ratios of 70 /xm / 24 /xm luminosities. Our analysis 
suggests that while star formation processes dominate the heating of the dust and PAHs, a heating 
process consistent with suprathermal electron heating from the hot gas, distinct from star formation, 
is heating the molecular gas and contributing to the heating of the ionized gas in the galaxies. The 
survival of PAHs and dust suggests that dusty gas is somehow shielded from significant interaction 
with the X-ray gas. 

Subject headings: galaxies: clusters: intracluster medium, galaxies: starburst, dust, infrared: galaxies. 



^ Perimeter Institute for Theoretical Physics, 31 Caroline St. N., Waterloo, Ontario, N2L 2Y5, Canada 
^ Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 



2 



Donahue et al. 



1. INTRODUCTION 

Infrared spectroscopy provides critical clues about the power sources of luminou s galaxies whose energy sources 
are s hielded from visual inspection by layers of dust a nd gas (iGenzel et alJ 119981 : iKennicutH 119981 : iLaurent et alJ 
I2000D . including star formation acti vity and AGN (e.g ., iRoussel et al.l 120011) . The Infrared Spectrograph (IRS) on 
board the Spitzer Space Telescope ()Houck et al.ll2004[ i. exploiting the sensitivity and spatial resolution of Spitzer, 
delivere d stunning infrared sp ectra from galaxies of many types. A key project, Spitzer Infrared Nearby Galaxy 
Survey (jKennicutt et al.ir2003| ). as a survey of nearby galaxies (D < 30 Mpc), was limited to ma .inly spirals and a few 
ellipticals. Studies targeting the brightest infrared galaxies in the sky (e.g., lArmus et al.ll200"9 ^ included mainly the 
most luminous IR galaxies (LIRGs) and nearby IR-bright, star-forming galaxies. Spitzer programs such as these have 
produced a treasure of infrared spectra of galaxies as well as improved and standardized techniques for measuring 
infrared features (e.g., iDale et all I2OO90 . These recent theoretical and observational efforts have identified useful 
infrared diagnostics, which now allow a physical interpretation of the spectra based on models (|Dale et al.l 120061 : 
IFarrah et al]l2007l : ISmith et al.ll2007bn . 

We explore here the infrared spectral signatures of Brightest Cluster Galaxies (BCGs). BCGs are not common 
even in large samples of galaxies, because massive clusters themselves are rare (of order a few in a box 100 Mpc on 
a side), so in general, to study an interestingly large and bright sample of BCGs, they must be specially targeted. 
Although most BCGs have red colors and are dust-free, suggesting little star formation, some 15-25 % show evidence 
of significant star formation (rates up to ~ 100 Mfr^ yr~^) in their UV and optical conti nua (e.g., I Johnstone et alJ 
1987; McNamara & O'ConncU 1989; Fa bianlfl99l ICrawford et al.|[l999l:lHicks fc Mushotik v 2005; Raffc rtv et al.ll2008l: 
Bildfell et al.1120081 : IHicks et al.ll20inHDonahue et al.ll20inll . 

Star-forming BCGs seem to be exclusively found in the centers of clusters whose hot intracluster medium (ICM) 

cores exhibit gas co oling times shorter than about 1 billion year s or low hot gas entro pies (K) where K — kTue '^^^ < 
30 keV cm-2, (e.g. iCavagnolo et al.l 120081 : iRaffertv et 311120081 : iHudson et all l201Cf) These clusters, once known as 
cooling flows, are called "cool-core" clusters. About half of the nearby X-ray luminous clusters fall into this category. 
This t rend for cooling flow clusters to host BCGs with powerful emission line nebulae was first found by iHu et ahl 
()1985[) . who noted a ~ 14 billion year threshold, limited by the far cruder X-ray data available at the time. Then, 
astronomers suspected that gas cooling from the hot phase was somehow related to these nebular emission line systems 
although the emission lines themselves were too bright to be generated directly by cooling gas. Almost two decades 
later, X-ray Multi-Mirror (XMM) spectroscopy failed to show [Fe XVII] and [O VII] emission lines, conclusively 
demonstra ting that very little X -ray emitting gas existed at temperatures 1/2-1/3 the temperature of most of the 
ICM (e.g., iPeterson et al.l l2003| ) . However, the disproof of the simplest massive cooling flow model did not explain 
why BCGs in these s ystems frequently exhibited properties indicative of activity: extende d emission line systems 
jHeckman et al."1989f), including vibrationally-excited molecular hydrogen at 1000 — 2 000 K l|Elston fc Malone"vlll994l : 
[Jaffc & Bremer 1997; Donahue ct al. 200(1), C O masses indicating c old R2 at - 100 K (lEdgdl200lD, UV excesses (most 
recently, O'Dea et al.ll2010t iHicks et al. 2010: iDonahue et al. 20101. radio sources (e.g.. Burns 1990: Cava gnolo et al.l 
I2OO8I IR emission from warm dust (Egami ct al.ll2006b : Donahue ct al. 2007b; O'Dea et al. 2008). It is important to 
note that the XMM spectral results ruled out the enormous X-ray cooling rates inferred from the simple cooling flow 
model 100s Mq yr~^), but do not provide limits near the cooling rates similar to the typical star formation rates 
estimated for these BCGs ('^ 1 — 10s Mq yr^^), with quantities an order of magnitude higher for the most extreme 
systems. Early Herschel results of a few classic examples of these active brightest cluster galaxies in cool-core clusters 
reveal far-IR spectra of similar sources that are consistent with the short wavelength Spitzer observations: strong 
peaks in broadband photometry from dust warmed by recently formed stars and powerful interstellar coolant lines of 
01 (63 ^m) and CII (153 ^m) (|Edge et al.ll2010ai rbl). and Spitzer IRS measurements of individual BCGs reveal that 
some, like Zwicky 3146 , have n ot only have powerful IR emission from warm dust but unusually luminous molecular 
hydrogen (Egami et al."2006a'), while others, like NGC4696, have luminous molecular hydrogen but only faint dust 
emission (Kaneda et al. 2008). 

These brightest cluster galaxies also pose a challenge to galaxy formation models. The so-called "over-cooling" 
problem in galaxy formation simulations creates massive galaxies that are bluer, even more luminou s, and with higher 
star formation rates than observed (see iBalogh et al.l 120011 : iCroton et al.l 120061 : 1 Bower et al.l I2006D . To remedy this 
situation, models must include AGN feedbac k in the for m of non-radiative energy in addition to stellar feedback to 
quench the formation of stars (e.g., ^S pi'ingc l et al.ll2005t ) and prevent the rapid cooling of hot intergalactic gas (e.g., 
iChurazov et~alll2001t iMcNamara fcNulsenI 120071 ) . Furthermore, the accretion of hot gas has been proposed as the 
dominant mode for forming the most massive (> 10^^'^ Mq) galaxies (e.g., iKeres et al.ll2005l ). 

Conveniently, BCGs provide a laboratory for this type of galaxy formation. Chandra observations clearly show AGN 
interactions with the hot ICM in the form of cavities in the hot atmospheres of cluste rs and galaxies. The mechanical 
energy associated with these cavities is is sufficient to offset cooling (Birzan et al. 2 0041 : iDunn fc Fabianl|2006l ), and 
to somehow regulate or quench star formation in most systems ( Raffertv et al.i 2006'). However, energetic feedback 
from AGN is apparently unable to offset cooling enti rely in all BCGs, and these systems that are struggling to offset 
rapid cooling are rich in cold gas and star formation (jO'Dea et al. 2008). The correlation between star formation and 
short cent ral cooling times in the hot gas show s that the gas fueling star formation may well have cooled from the 
hot ICM (|C avagnolo et al.l l2008t IRaffertv et ahl [2008) . Assessing and decoding the state of gas, dust, and stars in 
these galaxies, using Spitzer spectra, will yield clues about which physical processes are most relevant in determining 
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TABLE 1 
Observation Log. 



Cluster Rcdshift^ Spitzer IRS Mode Obs Date Duration (s) # of Slit Positions 

BCG ID SL LL SL LL 



2A0335+096 0.0347 20345 Staring 2006 - 09 - 16 10156 1698 12 12 

Abell 478 0.0860 20345 Staring 2006 - 03 - 17 10156 1698 12 12 

Abell 1068 0.1386 3384 Mapping 2005 -04- 21 2925 234 6 2 

Abell 1795 0.0633 3384 Mapping 2005 - 02 - 07 11701 1879 24 16 

Abell 1835 0.2520 3384 Mapping 2005 - 02 - 13 2925 234 6 2 

Abell 2597 0.0821 3384 Mapping 2005 - 06 - 30*^ 6826 704 14 6 

Hydra A 0.0549 3384 Staring 2005 - 12 - 14 7313 1509 12 12 

MS 0735.6+7421 0.216 20345 Staring 2006 - 04 - 25 21280 2768 4 4 

PKS 0745-19 0.1028 20345 Staring 2006 - 05 - 16 10640 1384 4 4 



^ BCG rodshift sources, fro m emission lines: 2A 0335+096 lIDonahue et al.ll20073). Hydra A ISmith et aU[200l). 
Abel l 1795 (CGCG 162-010) (Hill & OcEcrlc"199?), Aboil 2597 (PKS 2322-12) jVoit &: i3onahucjll997HColl ess eTaTl 
[200^ ) . seealsohttpj_/ /www. mso.anii.edii.au/2dFGRS/ Aboil 478 (NVSS 1041325+102754) I Zabludoff ct al. 199C|, PKS 
0745-19 I IHunstead"et al.iil978i), Abel l 1068 (FIRST J104044. 4+395712) lAllen et al..,1992,) . MS0735 (Z wCl 1370 or 
BCG 40 +74.12) llStocke et al.lll991l) . Abell 1835 (SDSS J140102. 07+025242. 5) (SDSS DR2 ; sec also lAllen et al.l 
IT99I) . 

^ Also 2005-07-05. 



TABLE 2 

Spitzer MIPS/IRAC Observation Log. 



Cluster 
BCG 


IRAC AOR 


Obs Date 




Duration 
(s) 


MIPS AOR'' 


Obs Date 




Duration'' 
(s) 


2A0335+096 


18646528 


2006 - 09 - 


26 


108 


18636544 


2007 - 02 - 27 




400, 300 


Abell 478 


11579904 


2005 - 09 - 


16 


1200 


14944256 (1), 14944512 (2 & 3) 


2006 - 02 - 22, 2006 - 03 


- 02 


550, 1500, 1200 


Abell 1068 


18650368 


2006 - 12 - 


27 


108 


18638336 


2006 - 12 - 08 




400. 320 


Abell 1795 










8788480 


2004 - 07 - 11 




36, 42 


Abell 1835 


4404480 


2004 - 01 - 


16 


3600 


4764160(1), 4744448(2 & 3) 


2004 - 02 - 20, 2005 - 06 


- 28 


1800, 600, 150 


Abell 2597 


13372160 


2005 - 11 - 


24 


3600 


13371904 


2005 - 06 - 18 




140, 150, 60 


Hydra A 


26923008 


2008 - 06 - 


09 


3600 


4707584 


2004 - 05 - 04 




140, 120, 180 


MS0735.6+7421 


7858688 


2003 - 11 - 


20 


500 










PKS 0745-19 


18667776 


2006 - 12 - 


27 


108 


18667520 


2006 - 12 - 08 




400, 300 



^ Channel (1) is the 24-/im channel, Channels (2 & 3) are the MIPS 70- and 160-/im channels respectively. 



Two durations indicate total exposure time in seconds for 24- and 70-/^m MIPS observations, respectively. Three durations indicate exposure times for 
24-, 70-, and 160-/^m MIPS observation sequences, respectively. 

a system's appearance and star formation rate during ac cretion of hot gas. 

We present here Spitzer Infrared Spectrograph (IRS) (jHouck et al.l 1200^ measurements of a sample of nine BCGs 
residing in cool-core clusters. We describe the measurement procedures, including the scaling we applied to the data 
to match aperture photometry, in §2. We present our results in §3. In §4, we compare the full spectra to a set 
of simulated time-averaged starburst spectral energy distributions (SEDs). In §5 we compare the emission line and 
polycyclic aromatic hydrocarbons (PAH) ratios and correlations seen in our sample to those seen in other types of 
galaxies. In §6 we discuss the correlation and lack of correlation in the various spectral components that suggest that 
at least two sources of heat must be considered in order to interpret the observations of these systems. We assume 
Ho = 70 km s~^ Mpc~^, and a flat, — 0.3 cosmology throughout. 

2. OBSERVATIONS AND DATA REDUCTIONS 
2.1. Observations 

The Spitzer IRS observations took place in 2005 and 2006 (see Table [T]), and the data were reprocessed in April 2009 
(vl8.7). Both short (SL) and long (LL) wavelength observations were obtained, at low spectral resolution, R ~ 60 — 130. 
With two spectral orders each, we obtained a total of four spectral modules. We have sparse spectral maps of nine 
cool-core BCGs, though we analyze only the central region here. We used IRSCLEAN vl.7 to apply the bad pixel 
mask supplied by the Spitzer pipeline and to find additional rogue pixels using a WCLEAN formula with an aggressive 
level of 0.5 (suitable for relativ ely faint targets such as these). The LL pixels are 5.1" across, while the SL pixels are 
1.8" across (|Houck et al.ll2004D . 

To cross-check our flux calib ration of the IRS spectroscopy, we also analyzed photometry data from the Spitzer 
Infrared Array Ca mera (IRAC) (|Fazio et al.ll2004[ ) and the from the Multiband Imaging Photometer for Spitzer (MIPS) 
(jRieke et al.ll200^ . We list the archival observations (known by their Astronomical Observation Requests or AORS) 
in Table [2j The aperture photometry is discussed in § 12.31 



2.2. Spectral Data Filtering and Extraction 



4 



Donahue et al. 



TABLE 3 

Parameters used in data reduction. 



Cluster Module FWHM FWHM Type Aperture center (J2000) Aperture size d'^ Factor*^ 







\ > 


(^Kpcj 






RA 


DEC 


\ ) 


V ) 




2A0335+096 


SL2 






Extended 


3 : 


38 : 


40.5 


9 : 


58 ; 


12 


10.8 X 10.8 


17 


1.16 




SLl 


8.0 


5.7 




3 : 


38 : 


40.6 


9 : 


58 ; 


11 


10.8 X 10.8 


17 


1.16 




LL2 


10.5 


7.5 




3 : 


38 : 


40.4 


9 : 


58 


9 


30.6 X 15.3 


-80 


1.01 




LLl 








3 : 


38 : 


40.5 


9 : 


58 ; 


10 


30.6 X 15.3 


-80 


1.01 


Abell 478 


SL2 






Extended 


4 : 


13 : 


25.4 


10 : 


27 


55 


7.2 X 7.2 


14 


2.26 




SLl 


7.0 


12.3 




4 : 


13 : 


25.4 


10 : 


27 


56 


7.2 X 7.2 


14 


2.26 




LL2 


8.4 


14.7 




4 : 


13 : 


25.3 


10 : 


27 


57 


15.3 X 10.2 


8 


1.71 




LLl 








4 : 


13 : 


25.3 


10 : 


27 


56 


15.3 X 10.2 


8 


1.71 


Abell 1068 


SL2 






Point 


10 


40 


: 44.5 


39 : 


57 


11 


12.6 X 5.4 


-40 


0.65 




SLl 


3.7 


10.6 




10 


40 


: 44.4 


39 : 


57 


11 


19.8 X 3.6 


-40 


1.00 




LL2 


9.0 


25.6 




10 


40 


: 44.7 


39 


57 


; 9 


35.7 X 10.2 


43 


1.00 




LLl 








10 


40 


: 44.6 


39 : 


57 


10 


56.1 X 10.2 


43 


1.00 


Abell 1795 


SL2 






Extended 


13 


48 


: 52.5 


26 : 


35 


33 


12.6 X 3.6 


-6 


1.87 




SLl 


8.1 


10.5 




13 


48 


: 52.5 


26 : 


35 


34 


12.6 X 3.6 


-6 


1.87 




LL2 


12.7 


16.5 




13 


48 


: 52.2 


26 : 


35 


35 


20.4 X 10.2 


78 


1.33 




LLl 








13 


48 


: 52.3 


26 : 


35 


36 


20.4 X 10.2 


78 


1.33 


Abell 1835 


SL2 






Point 


14 


: 1 


: 2.1 


2 : 


52 : 


41 


16.2 X 3.6 


-11 


1.20 




SLl 


3.8 


19.4 




14 


: 1 


: 2.1 


2 : 


52 : 


42 


16.2 X 3.6 


-11 


1.20 




LL2 


9.5 


49.2 




14 


: 1 


: 1.9 


2 : 


52 : 


40 


45.9 X 10.2 


72 


1.00 




LLl 








14 


: 1 


: 2.0 


2 : 


52 : 


39 


45.9 X 10.2 


72 


1.00 


Abell 2597 


SL2 






Extended 


23 


25 


: 19.9 


— 12 


: 7 


27 


9.0 X 3.6 


27 


3.42 




SLl 


6.4 


10.7 




23 


25 


: 19.8 


— 12 


: 7 


27 


9.0 X 3.6 


27 


3.42 




LL2 


15.0 


25.3 




23 


25 


: 19.7 


-12 


: 7 


26 


15.3 X 10.2 


21 


1.87 




LLl 








23 


25 


: 19.8 


-12 


: 7 


26 


10.2 X 10.2 


21 


1.87 


Hydra A 


SL2 






Point 


9 


18 


: 5.7 


-12 


: 5 


43 


14.4 X 3.6 


-20 


1.70 




SLl 


4.7 


5.3 




9 


18 


: 5.7 


-12 


: 5 


44 


18.0 X 3.6 


-20 


1.70 




LL2 


8.7 


9.7 




9 


18 


: 5.4 


-12 


: 5 


43 


35.7 X 10.2 


63 


1.70 




LLl 








9 


18 


: 5.5 


-12 


: 5 


44 


35.7 X 10.2 


63 


1.00 


MS 0735 


SL2 






Extended 


7 ; 


41 : 


44.6 


74 : 


14 


39 


5.4 X 3.6 


-7 


3.58 




SLl 


6.2 


27.5 




7 ; 


41 : 


44.6 


74 : 


14 


39 


3.6 X 3.6 


-7 


3.58 




LL2 


10.4 


46.0 




7 : 


41 : 


44.2 


74 : 


14 


39 


10.2 X 10.2 


-13 


1.72 




LLl 








7 ; 


41 : 


43.8 


74 : 


14 


39 


10.2 X 10.2 


-13 


1.72 


PKS 0745-19 


SL2 






Point 


7 ; 


47 : 


31.4 


-19 


: 17 


: 39 


12.6 X 3.6 


-11 


3.36 




SLl 


4.9 


10.3 




7 ; 


47 : 


31.4 


-19 


: 17 


: 36 


27.0 X 3.6 


-11 


1.29 




LL2 


9.4 


19.7 




7 : 


47 : 


31.3 


-19 


: 17 


: 37 


25.5 X 10.2 


72 


1.00 




LLl 








7 : 


47 : 


31.2 


-19 


: 17 


: 36 


25.5 X 10.2 


72 


1.00 



(a) Anglo of longer axis of aperture, in degrees east of north (CCW). (b) Factor applied to the extracted spectrum, to scale it 
up to the total light in a similar broadband aperture (only for extended-source spectra) and in a few cases to co-register modules 
and improve agreement with broadband measurements. 

This section discusses our choices regarding the IRS spectral extraction process, with particular atte ntion to the 
treatment of extended sources compared to point sources. We used the software package CUBISM vLlQ ([Smith et alj 
l2007ai) to combine the exposures into a data cube, or spectral map. To reduce noise near the ends of the slit, we 
trimmed the exposures in the cross-dispersion direction by 3 — 5%. We confirmed that the off-target, paired observations 
for each spectrum were indeed source-free. These blank sky spectra were used for background subtraction. 

To remove rogue pixels which were not caught by IRSCLEAN, we composed bad pixel lists using CUBISM's au- 
tobadpix algorithms, on both the global and individual record levels. At the global level, we flagged any pixel which 
deviated by more than 2.5(7 from the median level in at least 50% of its appearances in the cube. This method flags 
only a few pixels, but each of these pixels has a relatively large effect on the spectral map. At the record level, we 
conservatively flagged any pixel which is a 5a outlier in at least 75% of its appearances in the cube. This method flags 
more pixels, but only in individual exposures. We also manually removed obvious rogue pixels. 

The Spitzer pipeline is optimized for single- slit observations of point sources, including a correction for the light lost 
from the slit, which can be as much as 36% (jSmith et al.ll2007"al Figure 4). However, CUBISM includes an option to 
remove this slit-loss correction factor, in order to extract spectra of extended sources. To determine which targets to 
treat as point sources, we used CUBISM to create a map combining all SLl wavelengths, and averaged the two rows 
which covered the peak source emission. Using this light profile, we measured the FWHM along the slit (see Table 
[3]). We did the same for the LL2. These modules were selected because their signal-to-noise is highest. Spitzer's PSF 
has an average FWHM of 2.6" in the SL and 6.6" in the LLB Four of our nine galaxies have a FWHM in the SLl of 
< 5.5" and were considered to be point sources. The other five sources do not uniformly fill the Spitzer slit, but they 
are not well characterized as point sources. 

For the point sources, we extracted spectra using an aperture that just spans the slit (2 pixels) and a length 
that captures most of the light along the slit (see Table [3]). Sp itzer data are calibrated for point sources and for 
this kind of aperture. Single-pointing software such as SMART0 (jHigdon et al.|[2005 iLebouteiller et al.ll2010l) uses a 

^ http: //ssc ■ spitzer ■ caltech. edu/dataanalysistools/tools/c ubism/] 

http: //ssc . spitzer . caltech. edu/irs/irsinstrumenthandbook/T 
® http: //ssc . spitzer . caltech. edu/dataanalyslstools/tools/ contributed/irs/smartT] 
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similar aperture with its "tapered column" extraction, increasing the length of the aperture as the PSF broadens with 
wavelength, and further optimizing the extraction by weighting each pixel by its signal-to-noise ()Lebouteiller et ahl 
iOlO). In order to use the same software for all of our spectra, we used CUBISM to extract our point-source spectra, 
approximating the tapered column type of aperture. (For PKS 0745-19 SL2, we truncated the aperture to avoid a 
noisy region. For Abell 1068 SL2, CUBISM spreads the light from those two pixels across three rows.) To check 
our procedure, we also extracted the point-source spectra using SMART and found agreement to within about 10%, 
sometimes to within 2%. 

For the extended targets, we removed the pipeline slit-loss correction factor. Our apertures include much of the 
available light with good signal-to-noise. However, our sparse spectral maps do not cover the full extent of the source. 

After extraction, any noisy edges were trimmed from each order. In the observed frame, the four low-resolution IRS 
modules span the following wavelength ranges: SL2: 5.2 — 7.6 /im; SLl: 7.5 — 14.5 fim; LL2: 14.3 — 20.6 /im; LLl: 
20.5 — 37.5 fim. Spectra extracted from CUBISM are reported in units of MJy sr~^, so the size of the extraction 
aperture is used to convert spectra to units of flux density (mJy) . We then corrected to the rest frame for each target 
(see Table [Ij by dividing both the wavelength and fluxes by a factor of {1 + z). 

2.3. Aperture Photometry and Systematic Uncertainties 

The light collected by the narrow Spitzer slit and our sparse spectral maps represents only a portion of t he MIR light. 
Therefore, to obtain meaningful luminosities, we rely on IRAC and MIPS photometry (see Table|4]) (as in lEgami et al.l 
(|2006aD ). The broadband aperture diameters are approximately three times the source's FWHM measured using the 
IRS (or equal to the FWHM of Spitzer 's PSF for the 70 and 160 /J,m points) with subtraction of background computed 
from a larger annulus. We adjusted the aperture size to exclude unrelated foreground or background sources, and 
applied the suggested aperture corrections given by the MIPS Instrument Handbook. At 24 /im, the corrections are 
1.17 (for apertures of 26" and 30") and 1.13 (for apertures 50"). For 70 and 160 /im the corrections are 1.22 and 
1.752 respectively. The photometric uncertainty is 5% for IRAC, and 10, 20 and 20% for the MIPS 24, 70 and 160 /im 
points respectively. 

Because the spectra within the cited IRAC and MIPS circular apertures may differ from the spectra obtained from 
within our smaller, rectangular IRS apertures, our analysis and conclusions rely most heavily on relative quantities, 
i.e., ratios, rather than absolute quantities. All correction factors for each module are listed in Tabled The interested 
reader can recover the flux in the apertures listed in Table |3] by dividing the fluxes published here by this factor. 

In three cases (Hydra A, A1795, MS0735), we do not have complete IRAC and MIPS coverage in the IRS wavelength 
range, and we needed a robust, standalone scaling procedure. For this purpose, we developed a scaling procedure that 
did not rely on IRAC or MIPS photometry. We validated this procedure, with scaled IRS spectra for sources with 
IRAC and MIPS photometry, by comparing the resulting spectrophotometry to IRAC and MIPS photometry. We now 
describe our scaling procedure. 

For the sources we identified as extended, we fit Gaussian, azimuthally symmetric light profiles (and neglected the 
fainter extended haloes) to estimate how much of the source was included within a circle of the same area as our 
rectangular aperture. The spatial profile from the SLl map was used as the reference to determine the scale for the 
SLl and SL2 spectra for light from outside the rectangular aperture. The profile from the LL2 map was used for both 
LL orders. For point sources, this scaling has already been performed by the pipeline (See § 12. 2p . 

An additional scale factor is needed to match the orders of the IRS spectrum. Flux mismatches are expected between 
the spectral orders extracted with CUBISM, because it is impossible to perform an exact tapered column extraction 
using CUBISM. Even tapered column extractions with SMART sometimes show mismatches. In our sample, three 
targets had mismatches between the first and second orders (SL1/SL2 or LL1/LL2), and three had mismatches between 
the SL and LL. Some spectra are plagued by noise or decreased signal at the module interface (e.g. 2A0335 near 7.1 
/im). We selected the LL2 spectrum as the photometric reference point because of the good agreement between the LL 
spectrum and MIPS photometry in a simi lar aperture and because the LL is less vulnerable to slit loss due to pointing 
errors (because of its larger pixels) (Smith et al. 2007a). As mentioned above, the LL2 has superior signal-to-noise to 
the LLl. We used low-order polynomial continuum fits to match the modules when major features did not interfere. 

In three cases an additional overall factor was needed (Hydra A (1.70), A478 (1.55), MS0735 (1.12)), probably 
because of extended halo light contained within the broadband aperture but not represented by our Gaussian profile. 
This procedure worked well, because agreement between the IRS spectra and broadband photometry was relatively 
good, almost always within 10%, and usually within 5% (see Figures [T][2]) . Therefore, while the IRS spectrum appears 
to exceed the IRAC value at 8 /xm for Hydra A, the agreement between the integrated spectrum and the IRAC 
photometry point is actually excellent. Note that the comparison here is between the MIPS and IRAC photometry 
and the observer-frame IRS spectra integrated over appropriate bandpasses. These flux points should not be confused 
with the integrated IRS photometry in the rest-frame 24-/im MIPS bandpass reported in Table [5l 

For analysis and plots, we combine the uncertainty calculated by the SSC pipeline in quadrature with a 15% 
systematic uncertainty to all fluxes to account for the uncertainty in scaling the IRS photometry relative to Spitzer 
broadband photometry. The systematic uncertainty dominates in almost all measurements. An additional 5% absolute 
photometric uncertainty is applied when making comparisons with data from other telescopes 0. 

2.4. PAHFIT Spectral Decomposition 
® http: / /ssc. spitzer. caltech.edu/spitzcrmission/missionovcrvicw/spitzertelescopehandbook/ 
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TABLE 4 

Broadband Spitzer, Photometry 



Band 2A0335+096 Abell 478 Abell 1068 Abcll 1795 Abcll 1835 Abell 2597 Hydra A MS 0735 PKS 0745-19 



3.6 


9.2 


12 


4.2 


20 


2.1 


10 




2.5 


14 


4.5 


20 




1.4 


20 


3.0 


14 


4.5 


5.6 


12 


2.8 


20 


2.0 


10 




2.0 


14 


2.9 


20 


2.9 [14] 


1.0 


20 


2.1 


14 


5.8 


4.4 


12 


2.0 


20 


2.7 


10 




1.3 


14 


1.8 


20 


0.7 


20 


1.7 


14 


8.0 


3.5 


12 


1.7 


20 


7.5 


10 




4.5 


14 


1.9 


20 


4.1 [14] 


0.4 


20 


2.3 


14 


24.0 


2.4 


26 


1.6 


30 


74.8 


30 


1.8 [40] 


17.8 


30 


2.1 


50 


9.1 [30] 






10.2 


30 


70.0 


77.1 


70 


62.8 


70 


894.5 


70 


37.2 [70] 


175.0 


70 


89.0 


70 


155.2 [70] 






154.3 


70 


160.0 






56.4 


80 






317.0 


80 


42.0 


80 


181.8 [80] 











Note. — IRAC and MIPS waveband centers are in units of fim. Observer-frame fluxes are in units of mjy, and aperture diameters 
in arcseconds are given in brackets. Photometric uncertainties are ~ 5% for IRAC and, for MIPS, 10, 20, 20% for 24, 70, and 160 /im 
respectively. 



We used the spectral-decomposition pacliage PAHFIT vl.2 ([Smith et al.ll2007bl) to make empirical fits to the IRS 
spectra and to facilitate direct comparison with results from other workers using the same method. The short- 
wavelength PAHFIT results are plotted in Figures [T][2] and the full- wavelength results are shown in Figure |3l PAHFIT 
fits the following components: a starlight continuum, several thermal dust continuum components, broad PAH emission 
bands, narrow atomic and molecular emission lines, and broad silicate absorption bands (Figures [T][51) We customized 
the list of fitted emission features to a limited set, excluding those that were very weak. In the case of 2A0335, small 
parts of the spectrum (below 5.3^m and between 7.05 and 7.35/im) were excluded from the PAHFIT analysis because 
noise in those parts of the spectrum hindered a successful fit. 

3. RESULTS 

Our BCG galaxies exhibit a number of emission features from PAHs, ions, and H2 molecules. The H2 features are 
unusually prominent. The BCG spectra qualitatively fall into two general categories. Four galaxies (A1835, A1068, 
PKS0745, and Hydra A) exhibit the strongly rising IR contin uum at > 2 5 /im an d distinct PAH features characteristic 
of galaxies with strong signatures of star formation (e.g., Br andl et al. 2006; S mith et al . 2007b). In the remaining 
five cases, the 5-7/im continuum is dominated by cool stars, in contrast to the spectra of starbursts. 

Several spectral features, which do not correspond to a known emission feature, are artifacts of noise or the data 
reduction process. For example, in the spectrum of A1835, there is noise on the red shoulder of the 11.3 fj,m PAH 
band where the SL and LL modules do not perfectly align; this is also responsible for noise near 13 /im in A478. The 
junction between the LL2 and LLl accounts for some of the noise near 17 /im for A1835, and near 20 /im for 2A0335. 
A feature near 24 fim in the spectrum of AIO68 may be attributed to [Ne V]24.3 /im, as discussed in i; 15. 5[ but is 
probably spurious. An emission feature near 4.9 /im in A478, and possibly A1795, A2597, and PKS0745, might be 
ascribed to [Ar V]4.93 /zm, or to an unidentified PAH. The feature at 20.7 ^m in A478 is unidentified. The noise 
at 15.0 /zm in the spectra of A2597 and A1795, and 18.0 fim in Hydra A, appears to be spurious. Note that noise 
increases dramatically past about 33 /xm in the observed frame. 

The features we will examine most closely in this paper are the relatively bright forbidden emission lines of [Ne II] 
at 12.8 /im, [Ne HI] at 15.6 fim, and the PAH complexes at 7.7, 11.3, and 17 /im. We include an analysis of the 
correlation of the intensities of the brightest rotationally-excited molecular hydrogen transitions, S(2) and S(3), with 
those of other spectral features. A selection of line measurements and intrinsic luminosities including 1-a statistical 
errors from PAHFIT are presented in Table [5l The continuum fluxes and luminosities [vL^) are found from the 
feature- free continuum (the stellar blackbody, thermal dust components, and silicate absorption). The continuum 
measurements were determined by averaging across bandpasses 1 /im in width at 6 and 15 /im, and weighting by the 
24-/im JVIIPS response. All fluxes and luminosities in Table [S] are presented in the rest frame, at the rest wavelength. 
To recover observed fluxes, multiply by (1 -I- z). A more detailed analysis of the molecular hydrogen line ratios and 
excitation diagrams is deferred to a paper in preparation. 
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Fig. 1. — Detailed decompositions of nine cool-core galaxy clusters from 4.3 to 14 fim, utilizing PAHFIT l|Smith et al.|[2007al 1. Red 
lines represent thermal dust components; magenta, the stellar continuum. Their combination is a thick gray line. Broad PAH emission 
complexes are plotted in blue, and the unresolved emission lines arising from low-ionization or molecular hydrogen emission are plotted in 
violet and labeled at the top. The full spectral extraction is indicated by the green line, plotted over the rest-frame flux intensities and 
statistical uncertainties. In the two cases where the empirical PAHFIT detected silicate extinction, the extinction curve is represented with 
a dotted line using the axis at right; all components are diminished by the extinction. Appropriately transformed IRAC photometry is 
indicated by cyan squares (see Table [Sjl. 
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Fig. 2. — Figure [T] continued. 
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Fig. 3. — PAHFIT results over the full wavelength coverage. The continuum (stellar + thermal dust) is represented by a thick gray 
line, and the full spectral extraction by a green line. The rest-frame flux intensities are plotted with statistical uncertainties. Note the log 
scaling exaggerates the uncertainty of the faint, long-wavelength continua for several targets. Appropriately transformed IRAC and MIPS 
photometry points are indicated by cyan squares (see Table [sj. 
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TABLE 5 

Line and Continuum Measurements 



Line 



2A0335 A478 A1068 A1795 A1835 A2597 HydraA MS0735 PKS0745 



Ne II 12.8 iim Flux 17.7 
Ne II 12.8 fim Flux Error ± 0.3 

Ne II 12.8 fira Lum 0.0499 

Ne III 15.5 Atm Flux 11.6 
Ne III 15.5 Atm Flux Error ± 0.7 

Ne III 15.5 Atm Lum 0.0328 

PAH 7.7 Atm Flux ■ ■ ■ 

PAH 7.7 Flux Error 
PAH 7.7 Lum 

PAH 11.3 Atm Flux 31.2 

PAH 11.3 Flux Error ± 0.7 

PAH 11.3 Lum 0.0881 

PAH 17 Atm Flux 32 
PAH 17 Flux Error ± 10 

PAH 17 Lum 0.090 

H2 SI Flux 38.7 

H2 SI Flux Error ± 0.9 
H2 SI Lum 0.109 

H2 S2 Flux 13.7 

H2 S2 Flux Error ± 0.3 

H2 S2 Lum 0.0386 

H2 S3 Flux 40.8 

H2 S3 Flux Error ±0.4 
H2 S3 Lum 0.115 

Cont 24 Atm Flux (mjy) 2.66 

Cont 24 Flux Error ±0.06 
Cont 24 Lum 0.937 

Cont 15 Atm Flux (mJy) 1.45 

Cont 15 Flux Error ±0.1 
Cont 15 Lum 0.817 

Cont 6 A^m Flux (mJy) 2.60 

Cont 6 Flux Error ± 0.03 
Cont 6 Lum 3.66 



17.4 


47 


13.0 


26 


23.2 


24.6 


0.95 


38.9 


lb 0.4 


± 2 


± 0.3 


± 4. 


± 0.9 


± 0.6 


± 0.2 


± 0.3 


0.319 


2.4 


0.126 


5.0 


0.386 


0.176 


0.129 


1.05 


8.2 


28 


5.63 


16 


21 


22 


< 0.5 


14.3 


lb 0.6 


± 1 


± 0.7 


± 2. 


± 2. 


± 1. 




± 0.3 


0.15 


1.45 


0.054 


3.1 


0.36 


0.16 


< 0.07 


0.385 




ouu. 






42 


lOU 




1U 


± 6 


± 8 


± 9 


± 6 


± 10 


± 10 




± 4 


0.49 


15.4 


0.27 


103. 


0.71 


1.3 


< 0.6 


1.2 


17 


91 


22 


132. 


16 


83 


4.6 


44.9 


± 1 


± 2 


± 1 


± 2 


± 3 


± 3 


± 0.3 


± 0.8 


0.32 


4.7 


0.21 


25.5 


0.27 


0.60 


0.63 


1.21 


^14 


260. 




< 36 


< 48 


42 


< 89 


25 




± 80 








± 20 




± 10 


< 0.25 


13 


< 0.26 


< 7.0 


< 0.80 


0.30 


< 0.12 


0.68 


16.8 


< 20 


18.0 




49 


10. 


1.8 


23.6 


± 0.7 




± 1. 




± 2 


± 2 


± 0.2 


± 0.5 


0.310 


< 0.9 


0.174 


< 2.2 


0.82 


0.075 


0.25 


0.635 


7.0 


10.9 


6.4 


5.1 


14.5 


4.0 


1.1 


9.7 


± 0.4 


± 0.7 


± 0.3 


±2. 


± 0.8 


±1. 


± 0.2 


± 0.3 


0.13 


0.56 


0.062 


0.99 


0.242 


0.029 


0.15 


0.26 


27.8 


28 


22 


9.6 


56 


11 


3.73 


27.0 


± 0.5 


-I- 1 


± 2 


± 0.6 


-I- 1 


-I- 1 


± 0.09 


± 0.3 


0.511 


1.4 


0.21 


1.9 


0.93 


0.077 


0.509 


0.726 


0.931 


91.0 


2.12 


25.4 


2.57 


9.0 


< 0.04 


11.6 


± 0.05 


± 0.1 


± 0.07 


± 0.1 


± 0.1 


± 0.06 




± 0.03 


2.14 


582. 


2.56 


614. 


5.35 


8.1 


< 0.7 


38.9 


0.520 


28.4 


1.04 


6.13 


1.07 


4.49 


< 0.1 


3.13 


± 0.1 


± 0.2 


± 0.1 


± 0.3 


± 0.3 


± 0.2 




± 0.06 


1.91 


291. 


2.01 


237. 


3.57 


6.43 


< 2.73 


16.8 


1.14 


2.23 


1.46 


0.99 


1.0 


2.20 


0.36 


1.18 


± 0.05 


± 0.05 


± 0.04 


± 0.09 


± 0.1 


± 0.05 


± 0.02 


± 0.05 


10.5 


57.1 


7.04 


96. 


8.4 


7.8 


25. 


15.9 



Note. — Emission line flux values are in rest-frame units of 10 W m ^, while continuum ('Cont') fluxes are quoted in 
mJy. Luminosities are in units of 10^^ erg s~^. Fluxes, luminosities, and wavelengths are corrected to their rest-frame values 
(To recover the observed flux, multiply by (1 -|- z).) The la flux errors are statistical only. A 15% systematic error (0.06 dcx) 
is included in the plots and in the analysis, as discussed in ^ 12.31 Upper limits are 3(7. 

^ The flt to 2A0335's 7.7fim PAH complex failed because of decreased signal near the expected location of the feature. We 
were unable to set an upper limit in this case. 
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4. BEST FIT STARBURST AND OLD STELLAR POPULATION SED MODELS 

We fit the IRS spectra of the BCGs with an ensemble of simulated spectral energy distributions (SEDs). These 
models allow us to estimate the total infrared luminosity and associated star formation rates and to identify differences 
between the SED of a star-forming galaxy and the observe d SEDs. We use th e suite of starburst models, including 
a wavelength-dependent attenuation template, describ ed inlGroves et alj (I2008D together with an SED of a 10 billion 
year old stellar population, derived using Starburst99 (jLeitherer et al.lll999[ ). iGroves et al.l ()2008l ) model the starburst 
SED as the time-integrated sum of distinct HII regions and the the photodissociation regions (PDRs) surrounding 
them, over a range of cluster ages and cluster masses. These models were used to reproduce the SEDs of typical 
template starbursts such as Arp 220 and NGC 6240. The models span 5 different metallicities (Z = 0.05, 0.2, 0.4, 1.0, 
and 2.0 solar). Metallicity affects the prominence of the PAH features and the dust-to-gas ratio. The models sample 
6 compactness parameters C which characterize the intensity of the stellar radiation at the HII region/PDR interface. 
More compact HII regions result in hotter grains. C determines the location (the "temperature") of the dust peak and 
thus controls the mid-IR emission peak, a feature that is not well constrained by IRS spectral coverage. 

Groves et al. (2008) provide models for five gas pressures {P/k — W^, 10^, 10^, 10^, 10^ K cm~'^), spanning the range 
of lower pressures in star-forming galaxies to the higher pressures in ULIRGs. However, the infrared SED between 
5-25 /im is insensitive to the pressure, except, at a marginal level, the forbidden lines, and we checked that changing 
the gas pressure does not change the other parameters. Since the ICM pressure in cool-c qre clusters is^ 10^ — 10' K 
cm~^, the fits reported in Table [7] sample the range of 10^, 10^, 10* K cm~'^. Finally, the IGroves et al.l (|2008D models 
were computed for unobscured HII regions and HII regions with PDRs. We opted here to fit only the PDR models, 
since the models with only HII regions had no PAHs. Models are normalized in luminosity to a star formation rate of 
1 M0 yr~^ sustained over the 10 Myr lifetime of the HII region. The spectra and the best-fit components are shown 
in Figures |4][5l 

We compared the observed spectra with the full suite of theoretical spectra and quantitatively identified the spectrum 
that best fit the data, in L^, units, by minimizing the sum of chi-squared. In addition to the quantized parameters 
described above, we fit the normalization for each starburst SED {Mq yr"'^), the attenuation {Ay), and the normal- 
ization of the old stellar SED (Mq). The Groves et al. (2008) models include forbidden lines and PAH emission, but 
not molecular emission, so we added a two-component local thermodynamic equilibrium (LTE) H2 spectrum, each 
component specified by T (K), and column d e nsity iV (10^* cm~^). We used molecular properties and A-values from 
iBlack &: Dalgarnol (|1976[ ): iHuber fc Herzbergl (|1979| ). The widths of the H2 lines were fixed to 0.05 /zm. 

In this comparison, the compactness parameter (C) seemed to divide BCGs with strong rising IR continuum from 
those with relatively fiat IR continuum in the IRS spectral range. The fit quality was sensitive to C, which mainly 
affects the steepness of the long-wavelength IR continuum. The BCGs with the best-fit C — 4 are 2A0335, A1795, 
A2597, A478, and MS0735. The first four of these galaxies have the lowest four 24/70 iim MIPS photometry ratios 
in the sample (= 0.03 ± 0.01); MS0735 lacks MIPS photometry. The signal to noise in the IRS spectrum of MS0735 
at long wavelengths is so low that there are no strong preferences for any starburst SED over another. The other 4 
BCGs in the sample have higher 24/70 fj,m ratios (= 0.08 ± 0.02). This trend in fiux ratios is consistent with the 
compactness parameter governing the peak in the dust spectrum, in the sense that more compact HII regions have 
PDRs with hotter dust. 

The metallicity also affects the SED in this region, particularly the relative strength of the PAH feature to the 
infrared continuum. The preferred starburst SED metallicities for the BCGs were at the high end {Z/Zq = 1 or 2). 
For all of the BCGs, the fit quality was very similar for Z/Zq = 1 or 2, except for that of A1068, which preferred 
Z/Zq = 0.5 or 1. That result is consistent with the fact that the BCG in A1068 has the lowest PAH to mid-IR 
continuum ratio of the sample (see § 15. ip . We note these high metallicities are consistent with the metallicity of 
the ICM in the centers of these clusters. As expected, pressure did not affect the fit quality in any significant way. 
The spectral features, including the 7.7 and 11.3 PAH features, were often well-matched by models including the 
PDRs at rest wavelengths shortward of 10-11 fim. The two-temperature molecular hydrogen template matched the 
H2 spectrum very well in most cases. 

To check th is process, we estima ted the total infrared luminosity from the expression Ltir = 1.559^24 + 2.1156^70, 
adapted from lDale fc Heloul ()2002[ ). L24 is i/L^, at 24 /im and L70 is vL^, at 70 fim. To adapt their equation 4 to this 
expression, we assume that vLi, at 70 /im ~ vL^, at 160 /im. Some assumption was necessary since the 160 /im flux 
is only available for 4 objects. This approximation may underestimate the TIR luminosity by about 30% compared 
to estimates using the 160 micron luminosities, but these 160 /im luminosity estimates are subject to large systematic 
uncertainties, and the detected sources may not be representative of the sources lacking photometry. (For MS0735-I-74, 
we have no MIPS photometry and only a 3-(t flux limit at 24 /xm from the IRS spectrum. We do not include this source 
in plots of TIR.) We compared Ltir e stimated f r om th is adaptation of Dale & Helou (2002) to Ltir derived using 
the calibration of L24 in Equation 2 of iWu et al.l (|2010[ ). To get rest-frame 70/24 /im ratios used in this calibration, 
we convert the observed MIPS luminosity 70 to 24 /xm ratios to rest-frame ratios for 8 of the 9 galaxies in the sample 
with k-correction factors (^79/24) based on the best-fit SEDs (< 3 — 25%; Table[6l). The derived total IR luminosities 
were consistent to better than 20% (Table |6]). The relationship between either estimate of Ltir. and the SED- inferred 
SFR is consistent with the Kennicutt (1998) SFR relation for starbursts (Figure [6^). This consistency indicates the 
24 and 70 /im MIPS data, not included in the fits, are consistent with the star formation rates inferred from the IRS 
data alone. 

In Figure [Hbj we plot the star formation rates based on the SED fits versus the lCalzetti eFaLl (|20Tot) mean relation 



12 



Donahue et al. 




10 20 30 

Wavelength (microns) 

Fig. 4. — Rest waveleng th 5-35 fim spectra, observed fluxes. (To convert to rest-frame flux, divide by (l-|-z).) Red line: Best flt constant 
star formation model from I Groves et aU 1)20081 ') . Blue line: Best fit old stellar population. Green line: Molecular hydrogen, two-temperature 
LTE model. 

for SFRs derived from the observed MIP S 70 /xm lu minosities {I'L^) (SFR(M0) = 0.059(L7o/10''2) ergs^i). The dotted 
hue shows the same relation, also from iCalzetti et al. (20l0), if a different Ha-based SFIl from Kennicutt (1998) is 
used for calibration. Figure |6] demonstrates reasonable consistency between SFRs estimated based on SED fits to the 
IRS spectra and SFRs derived from MIPS photometry. One caveat to this comparison is that the IR continuum of 
galaxies with a significant old stellar population may have contributions from dust heated by these cool stars. The 
observed 70-/im luminosities of the BCGs are somewhat higher than what is predicted by the best-fit SEDs, by a 
factor of 1.3-2.4, which, if taken literally may indicate that dust heated by evolved stars may be contributing between 
30-60% of the 70 /xm luminosity. T he 70 /im luminosities of these BCGs are in the valid domain for applying the 
relationships in ICalzetti et al.l ()2010t ): the 160 /im luminosities are likely to be even more contaminated from cool dust 
emission unrelated to star formation. 

It is interesting to note that the SFR estimates based on 24-/im luminosities, from Calzetti et al. (2010), are 
systematically lower compared to SFRs based on these SED fits or the 70-/xm luminosities for the 4 galaxies with high 
70/24 ^m flux ratios (A478, 2A0335, A2597, and A1795) and low best-fit C = 4, as well as MS0735 (which has no 
MIPS photometry). In Figure [7l we plot the rest-frame 70/24 /im luminosity ratio of the BCGs together with the 
local sample studied by Calzetti et al. (2010; Figure 17). The rest-frame 24 /xm luminosity is calculated from the 
IRS spectra, but it is consistent with MIPS photometry. The rest-frame ratios are based on MIPS photometry alone, 
with small k-corrections (< 10% for z < 0.1, and 10 — 25% for z > 0.1), based on the best-fit starburst SEDs. The 
70/24 /im ratio varies by about factor of 10 for a given 24 /tm luminosity in the full sample of Calzetti galaxies, and 
some of this scatter is due to metallicity, with the highest metallicity systems having the highest 70 / 24 //m luminosity 
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Fig. 5. — Figure |4] continued. 
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Fig. 6. — A comp arison of infrared Brightest Cluster Galaxy (BCG) star formation rates (SFRs) from the best-fit starburst models of 
I Groves et al.l I I200SI) and total infrared luminosity {Lxir) and the SFRs inferred from Lxm- On left, we show these inferred SFRs are 
highly correlated with the total infrared luminosity (Lxir DH; Table[6ll. For comparison, the solid line shows the Kennicutt (1998) relation 
for starbursts ( SFR(Mq yr"-*^) = 0.045(Ltir/10^^) erg s"-"^ ), along with ±30% typical calibration scatter (dotted lines) . On right, the 
same SFRs are plotted against the SFRs based on 70 /im luminosities, using the mean relation in [Calzetti et al. 1 12010) (SFR(M0) = 
0.05 9(L7o/10*^) ergs~^). The solid line shows the line of equality; the dotted line shows how the predicted SFR from L70 would change, 
from lCalzetti et aLI 1 12010) . if a different Ha-based SFR from Kennicutt (1998) is used for calibration. 
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TABLE 6 

Total Infrared Luminosity Estimates 



Name L24 Ltir (DH) Ltir (W) ^70/24 logL7o/L24 

10*2 erg 10*2 erg s'^ 10*^ erg 8"^ (rest) 



2A0335 


0.94 


21.2 


25.8 


0.97 


1.07 


A478 


2.14 


108. 


76.2 


0.92 


1.09 


A1068 


582. 


5060. 


4120 


0.73 


0.48 


A1795 


2.56 


36.6 


40.9 


0.94 


0.82 


A1835 


614. 


4028. 


4210. 


0.83 


0.45 


A2597 


5.35 


143. 


177. 


0.92 


1.13 


Hydra A 


7.32 


112. 


94.3 


0.94 


0.74 


PKS0745 


38.9 


437. 


389. 


0.84 


0.64 



Note. — All luminosities are in units 10^^ erg s ^ . L24 is rest-frame i^Lj^ at 24 /^m 
from the IRS spectr a. Ltir (DH) is the TIR derived from a relation adapted from 
pale &: Heloulf20Q^ . Ltir (W) is based on L24 and the observed ratio L70/L24 
luminosities based on MIPS photometry, from Equation 2 in W u ct a l. lj2010). The 
k-correction for the 70/24 /j-m luminosity ratio is based on rest- and observer-frame 
MIPS response functions convolved with the best-fit SED shapes. The last column is 
the logiO ratio of the k-corrected, rest-frame L70/L24 luminosities. 
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Fig. 7. — The L70/L24 ratio as a function of the 24 fim luminosity L24.The open black triangles are higher-metaUicity star forming 
galaxies; the closed blue triangles are lower-metallicity star-forming galaxies; and the open red squares are luminous infrared galaxies 
(LIRGs), all from Calzetti et al. (2010). The green symbols with error bars are the rest-frame 24 fira luminosities and rest-frame (k- 
corrected) 70 to 24 /im luminosity ratios for the 8 BCGs in our sample with MIPS photometry. 



ratios. The 70/24 /im ratios of BCGs are similar to, but somewhat higher than, those of high-metallicity star-forming 
galaxies. The two BCGs most like starbursts, A1835 and A1068, have ratios similar to the luminous infrared galaxies 
(LIRGs) of the Calzetti et al. (2010) sample. 

We briefly discuss the results of the comparison to theoretical SEDs of star-forming galaxies for each source, ordered 
approximately by decreasing mid-IR luminosity. Note that all fluxes and luminosities in this paper are derived from 
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TABLE 7 

Best Fit Starburst PDR + HII Model Parameters and Stellar Masses 



Name 




C 


Log P/k 


HII 


Av 


Old Star 


SFR 


Tl 


Nl (Ha) 


T2 


N2 (Ha) 








(K cm~^) 


or PDR 




igll M0 




(K) 


ClOl* rm-^) 


(K) 




2A0335 


1 


4 


8 


PDR 





1.5 


0.7 


350 (5) 


52 (0.2) 


1200 (30) 


0.8 (0.1) 


A478 


2 


4 


8 


PDR 





4.4 


2.7 


460 (10) 


13 (0.7) 


1920 (200) 


0.12 (0.02) 


A1068 


1 


6.5 


8 


PDR 


2.5 


8.6 


100 


320 (4) 


98 (7) 


2000* 


0.120 (0.007) 


A1795 


2 


4 


8 


PDR 





2.6 


2.3 


260 (13) 


48 (0.6) 


760 (40) 


1.7 (0.2) 


A1835 


2 


5.5 


6 


PDR 





26 


270 


530 (40) 


4.7 (1.5) 


1300 (750) 


0.03 (0.02) 


A2597 


2 


4 


8 


PDR 





4.1 


5.4 


240 (4) 


160 (8) 


810 (17) 


4.9 (0.3) 


HydraA 


2 


5 


7 


PDR 





3.3 


4.3 


380 (25) 


15 (0.3) 






MS0735 


2 


4 


6 


PDR 





9.7 


0.3 (0.12) 


300 (240) 


4.4 (7) 


730 (70) 


0.5 (0.2) 


PKS0745 


2 


6.0 


8 


PDR 





5 


11 


390 (20) 


26 (1) 


1200 (65) 


0.48 (0.05) 


Note. — 


*Valuo 


pegged 


at extreme 


temperature 


. Uncertainties are 


quoted in p 


arcnthcscs. 


Component normalizations for 


the old stellar 



population SED and the starburst SED had statistical uncertainties of less than 1-3%, with the exception of MS0735. However, the photometric 
calibration and scaling uncertainties were ~ 15 — 20%, so at minimum, uncertainties at that level apply to the stellar masses and SFRs estimated 
here. 

PAHFIT, not these fits. 

Abell 1835. The IRS spectrum of Abell 1835 shows a strong red continuum and PAH features whose shapes and 
intensities are well-fit by the model in Table [71 The most prominent residual is the overestimate of the flux of the 
[S III] 18.7 fim hue. 

Abell 1068. The IRS spectrum of Abell 1068 is similar to that of Abell 1835 in that it has a strong red continuum, 
but its emission features are weaker. All features are well matched except for a 24 /im feature that is likely spurious but 
could be associated with [Ne V]24.3 fim. (The deficiency in the model spectrum between 13-18 /xm is a shortcoming 
of the iGroves et alJ ()2008l ) models compared to many starburst spectra.) 

PKS0745-19. The IRS spectrum of PKS0745-19 also has a strong red continuum but fainter PAHs than Abell 
1835. The observed long wavelength slope is flatter than the model. The observed blend of [Fe II] + [O IV] at 26 /im 
is bright. For this source, the [Ne II] 12.8 /im and [Ne III] 15.6 /im features are underpredicted by our simulated SED. 
The H2 S(l) transition is quite bright. 

Hydra A. The fit to the IRS spectrum of Hydra A overpredicts the [S III] 33 /im line somewhat, while the [S III] 
18.7 /im line is underpredicted. The ratio of this line pair is set by gas density. Shortward of 20 /im, the spectrum 
resembles the best fit starburst model fairly well, except for the previously noted deficiency in the model spectra 
between 10-19 /im. As in other systems, the [Fell] -I- [OIV] blend is not modeled. 

Abell 2597. The IRS spectrum of Abell 2597 is fairly flat. Most features seem to be represented. There are small 
excesses of [Ne II] and [Ne III] emission visible in the residuals. 

Abell 1795. The IRS spectrum of Abell 1795 is fairly well represented by the starburst models. 

Abell 478. The IRS spectrum of Abell 478 is faint. The data exhibit stronger lines of [Ne II] and [Ne III] and a 
flatter long-wave spectral slope. 

2A0335+096. The IRS spectrum of 2A0335+096 has a flat red continuum, and molecular hydrogen and neon lines 
are well flt by the model. The most prominent feature in the residual spectrum is the strong [Si II] 34.8 /im line. 
2A0335 is the lowest redshift source, so we have no other sources for context, and it docs sit near the edge of the 
spectrum where noise spikes are not uncommon. The [Ne III] /[No II] ratio is higher than predicted in the best fit 
model. 

MS0735+74. The IRS spectrum of MS0735+74 has a fiat red continuum with very weak emission line features, 
and to the limits of the data, fit but not tightly constrained in this exercise. 

In summary, the starburst models of time-averaged HII regions and PDRs, based on fits to the IRS data, do a 
surprisingly good job at qualitatively representing the continuum, PAH, and nebular features of the IRS spectra of 
BCGs, but far from perfectly. On the other hand, the star formation rates derived are consistent with estimates based 
on the 70/im continuum or PAH features. We will discuss this further in 15.31 where we show that the H2 and [Ne II] 
luminosities are significantly higher than those of star-forming galaxies with similar infrared luminosities. While the 
PAHs and the IR continuum are usually well represented, the models do not match the nebular [Ne II] (and [Ne HI]) 
emission relative to the continuum, and the slope of the continuum through the longest wavelengths of the IRS spectra 
is not consistently fit for spectra with faint IR continuua. We note that the metric for a best-fit for a given starburst 
model is dominated by the continuum since most of the points arc continuum-dominated. 

5. DISCUSSION 

In typical star-forming galaxies, the luminosities of dust, PAHs, [Ne II] lines, and even the rotationally-excited 
molecular hydrogen lines are linearly correlated with each other and with the star formation rate (SFR). There is a 
significant uncertainty in the SFR inferred for any individual galaxy, a factor of 5-10, because of the dispersion, but 
these quantities are highly correlated in star-forming galaxies. Therefore we compare the correlations we see for the 
BCGs in our sample with those of star-forming galaxies. We will show here, based on the correlations and ratios that 
we observe for BCGs, that the infrared continuum and PAH features are consistent with being powered primarily by 
star formation in BCGs. In contrast, the emission lines from rotational transitions of hydrogen are uncorrelated with 
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Fig. 8. — PAH and infrared continuum properties compared. The 24 /im continuum luminosity (L24) is strongly correlated with the sum 
of the luminosities of the PAH complexes at 11.3 fim and 7.7 fim. The slope of the best-fit power law (plotted as a solid line) is 0.90 it 0.03. 
For comparison, a dotted line of slope unity is shown. On the right, we show that the ratio of PAH (11.3 iim) to Ltir (W; Table |6)| is 
~ 0.0039 ± 0.0020. The solid line is the mean ~ 0.00661q;oo42 from lWu et~an ||20TOD : the dashed line is the lower limit of the range of their 
best-fit to the mean for starburst galaxies. Given the uncertainties in converting from L70 and L24 to Lt/jJi this plot shows that these 
galaxies have PAH/IR luminosity ratios that are only somewhat lower than normal star-forming galaxies, with the exception of A1068 and 
A2597. 

the dust and PAH features, and are primarily powered by a second process. The forbidden lines of neon are correlated 
with the IR emission, but not linearly. This pattern is consistent with these lines being powered by star formation 
and a second process that does not provide much heat to the PAHs and dust but is very effective at producing H2 
emission. This second heating mechanism is consistent with heating by a population of suprathermal electrons, either 
from the hot gas or perhaps associated with the radio source. 

All linear correlation coefficients (r) in the analysis below are based on the measurements in logarithmic quantities 
unless stated otherwise. For reference, a correlation of r > 0.66(r > 0.86) might be considered significant at the 2a 
(3cr) level for iV = 9 points (e.g. Bevington 1969; Bevington & Robinson 2003). In general, if two quantities are 
correlated linearly (i.e., the log-log slope is unity), they may have common origins. But if they are correlated, but not 
linearly, there may be something more interesting going on. For that reason, we plot dotted lines with unity slope in 
our correlation graphs. 

5.1. Dust and PAH Luminosity Correlations 

The smallest dust grains are the polycyclic aromatic hydrocarbons (PA Hs), composed of only a hundred atoms or so . 
These structures generate e mission from C-H or C-C-C bending modes ( Leger fc Pugetlll984l: iBoulan ger et al.lll998l : 
I Van Kerckhoven et al.ll200"ol ). excited by the absorption of UV photons ( AUamandola et al.l ll985: Slo an et al.lll999D . 
Such photons can heat these tiny grains stochastically, causing them to suddenly increase in temperature then cool 
(e.g.,|Li & Drainc 2001). PAH features at 3.3, 6.2, 7.7, 8.6, and 11.3 /im in spectra, are th ought to be an excellent tracer 
of B stars, or of relatively recent star formation (jPeeters et al.ll2004t iBrandl et al.ll2006l : fPorster Schreiber et al.l [20041 . 
The continuum I R luminosity, wh ich traces star formation, and PAH luminosity is strongly c orrelated in norm a l star - 
forming galaxies ()Wu et al.ll2010D . Studies of low n ietallicity star-forming dwa rf galaxies by [Rosenberg et all (|2008[ ) 
and of star forming regions in irregular galaxies by iHunter fc KaufmanI ()2007D show the PAH emission decreases as 
metallicity decreases, so metallicity is one factor that can lead to scatter in the correlation between PAH emission and 
IR luminosity from dust. 

The brightest features from the polyaromatic hydrocarbons (PAHs) are the complexes at 11.3 fim and 7.7 fim. The 
sum of those lines in our BCG sample is strongly correlated with the 24 fim continuum in both flux and luminosity 
(Figure E]). The correlation coefficient r = 0.93 for both. Excluding MS0735 does not affect the correlation. The 



relationship is very close to linear: Lii.3+7.7 cx L24^°^° °'^, where L24 



vLi, at 24 fim rest frame. (L11.3 oc L'. 



0.96±0.05 
TIB. 



The ratio of the 11.3 fim luminosity to Ltir is 0.0039 ± 0.0020 (omitting MS0735 from the sample for lack of 
MIPS data), plotted in Figure [8b, to compare to the mean of ~ 0.0066±i^;[^^4 ^L°]^°^;^^°^ from 123 starburst dominated 
galaxies from a 24- /^m flux-limited sample of 330 galaxies in lWu et"aLl ()2010D The best-fit p ower law relating Ln.a^m 
is a = 1.05 ± 0.05, also similar to that seen for starburst galaxies (jWu et al.ll201 0V Given the uncertainties 
in converting from Ljq and L24 to Ltir, this comparison shows that these galaxies have PAH/IR luminosity ratios 
that are only somewhat lower than normal star-forming galaxies, with the exception of AIO68 and possibly A2597. 
The PAHs in Abell 1835 are about 4 times brighter compared to L24 than the nearly equally IR luminous Abell 1068, 
so we detect significant intrinsic scatter in this ratio. There is no correlation in the ratio of PAH/IR luminosities to 
IR luminosity (see Figure [8)d). 

The 11.3 fim PAH luminosity is highly correlated with the 24 fim and the TIR luminosity, which makes sense if the 
dust and PAHs are heated by the same process. The 15 and 24 fim continuum luminosities are strongly correlated 
with each other (r = 0.96) and nearly linearly correlated (L24 oc ii^042±0.006-j^ expected since both quantities are 
usually produced by dust grains. 

In contrast, 24 fim continuum and 11.3 fim PAH luminosities are not correlated with 6 fim continuum luminosity. 
In fact, when MS0735 is excluded, there is no correlation between 24 fim and 6 fim luminosities, r = 0.17. Similarly 
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Fig. 9. — A plot of the PAH 7.7 to 11.3 /xm ratio vs. the ratio of the 24/6 /xm continuum luminosities, The PAH ratio indicates the 
ionization level of PAHs, and this ratio is higher in systems in which the 24 ^m luminosity from dust exceeds the 6 /^m stellar continuum. 
This trend is consistent with the PAHs in the BCG systems most like starbursts experiencing a harder incident radiation spectrum. 

excluding MS0735, there is no correlation between the 6 /im and PAH flux at 11.3 /ini (r = 0.23) or [Ne II] (r = 0.004). 
(Including MS0735 in the tests increases the correlations to ~ 0.6, under the 2a threshold, but because the computed 
significance relies on the inclusion of a single source, it must be considered spurious.) The 6 /im light in these BCGs is 
produced primarily by old (cool) stars, and therefore is a metric for the stellar mass. The lack of correlation between 
dust and stellar continua luminosities suggests that dust, PAH, and gas heating is not determined by cool stars. The 
systems where the dust luminosity well exceeds the 6 /Ltm luminosity from stars, A1068 and A1835, exhibit higher 
PAH 7.7 to 11.3 /im ratios, consistent with the hypothesis that these are like starburst galaxies with levels of PAH 
ionization similar to those seen in starbursts (Figure ^ . We will consider these ratios more fully in § 15.41 

For the rest of this discussion, we will assume that the long wavelength IR continuum is powered primarily by 
obscured recent star formation. However, even though dust heating by evolved stars does not seem to dominate 
these systems at 24-70 /um, evolved stars may be the dominant source of heat for cooler dust emitting at longer 
wavelengths, and this dust therefore could contribute to emission at shorter wavelengths (see § Furthermore, 
processes of interest such as suprathermal electron heating and weak AGN may also supply energy to these systems. 
Since the observed global quantities are galaxy-wide averages, they suffer from the same interpretation ambiguitiy as 
high redshift, unresolved sources. Even if we could interpret these spectra in the context of star formation alone, it 
is impossible to unambiguously distinguish between a star formation episode of a single age and stellar mass and a 
time-averaged star formation history of "constant" star formation. Physically, signatures of star formation in normal 
galaxies include cold, dusty molecular gas, excess UV continuum, Ha, PAH emission, and infrared dust emission. We 
will discuss these data in a framework where obscured star formation is tracked by the infrared and PAH emission. 
However, we will show that star formation alone is inadequate to explain the full set of infrared spectral features in 
these systems. 



5.2. Forbidden Neon Line Correlations 

The luminosities of forbidden lines of neon, which are channels for radiative cooling, are sensitive to the thermal 
energy input into the ionized gas . They therefore h ave also been shown to be good tracers of star formation rates 
in normal star-forming galaxies. iHo fc Ketol (|2007| ) showed that the sum of the fine structure lines of Ne II (12.8 
/im) and Ne HI (15.6 /im) correlates strongly with IR luminosity in normal star-forming galaxies over 5 orders of 
magnitude in luminosity. The sums of [Ne II] and [Ne HI] luminosities in our BCG sample also correlate strongly 
with L24 (r = 0.95, fiuxes correlate with r = 0.90). [Ne II] alone is just as correlated (r = 0.94, fluxes correlate with 
r = 0.91). The relationship, however, deviates even more from linearity than the PAH-IR luminosity relationship, with 
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Fig. 10. — [Ne II] and infrared continuum properties compared. On the left, [Ne II] luminosity is strongly correlated with L24 luminosity, 
but the best-fit power-law (solid) is flatter than linear, with a index of 0.59 ±0.03. A dotted line with slope unity is plotted for comparison. 
On the right, the ratio of [Ne II] luminosity to Ltir {Ltir (W); Table |6ll continuum luminosity decreases with increasing Ltir- The ratio 
for typical star-forming galaxies, log([Nc IIJ/Ltir) = —3.44 ± 0.56, from Ho & Koto (2007) is plotted. The dotted line is the approximate 
upper bound of the intrinsic scatter. 

L{[Ne II]) scaling as iP^^^^^-^^ (or L^J^^^" ""*). We suspect that while the dust and the PAHs are heated primarily 
by star formation, this lack of linearity in the [Ne II] -IR correlation suggests that star formation may not be the sole 
process producing [Ne II] emission. 

The ratio of [Ne II] to total infrared lum inosity {Ltir (W) ; Table H]) decreases somewhat with increasing IR lumi- 
nosity (Figure 

Wi- The IHo fc Ketol (l2Q07l) mean relationship between [Ne II] and Ltir for star-forming galaxies is 
log([Ne II] /Ltir) = -3.44 ± 0.56, nearly independent of Ltir. Excluding MS0735, the [Ne II] luminosity in BCGs is 
about 1.6 — 12 times higher than the mean [Ne II] luminosities of normal star- forming galaxies of similar infrared lumi- 
nosities. The largest differences are found for the BCGs with lower IR luminosities (< 11^° L©)- MS0735-I-74, for which 
there is only an upper limit continuum estimate, is particularly bright in [Ne II] compared to its i nfrared lurn i nosity 
(> 0.009), > 25x the mean. The observed scatter of this ratio for normal star- forming galaxies in lHo fc Ketol (|2007l ) 
is large, ±0.5 dex; nevertheless, the BCG ratios sit consistently on the high side of the scatter for normal star-forming 
galaxies, indicating that another process beyond star formation is also contributing to the heating of the ionized gas, 
particularly in the low-luminosity systems. In summary, the [Ne II] luminosities seen in the low IR-luminosity BCGs 
exceed what would be expected from a star-forming galaxy with the same IR luminosity, but the two quantities are 
strongly correlated. 

The [Ne III]/[Ne II] ratio is not at all correlated with the mid-IR luminosity (Figure fTTj) . It is possible that the [Ne 
III]/[Ne II] ratio indicates an approximate starburst age, with BCGs having the highest [Ne I II]/[Ne II] also having 
the youngest starburst populations, but not necessarily the largest numbers of young stars (|Thornley et al.l 120001 : 
IRigbv fc Rieki[200l ISd^ders et al.l[2007l) . 

5.3. Molecular Hydrogen Correlations 

Extremely luminous pure rotational II2 lines, usually S(3) 9.67 /xm, S(2) 12.28 /im, S(l) 17.04 /im, are detected in 
all nine galaxies. S(0) 28.22 /im was not detected in any of these sources. Rotational transitions from S(5) - S(7) (5.51 
^m) are seen in a majority of these spectra. While the lumi nosities of rotation ally-excited molecular hydrogen lines 
are correlated with IR luminosities of star- forming galaxies (jTrever et al.ll2010l) . that is certainly not the case with 
our BCG sample. The line luminosities from rotational molecular hydrogen transitions from these BCGs are much 
greater than expected from the level of star formation heating the warm dust. II2 emission is also uncorrelated with 
the continuum at 15 or 24 /zm, r = 0.4 — 0.5. 

Rotational emission fr om molecular hydro gen is commonly detected in ULIRGs ( Higdon et al.l [20061 ) and in star- 
forming normal galaxies (jRoussel et al.l[2007[ ). In such galaxies, the luminosities of these lines are only about 4 x 10"'* 
of the total infrared power between 8-1000 /xm. However, the ratio of II2 luminosity to L24 for the BCGs in this sample, 
ranging from 0.004-0.3, is about 5-100 times more than one would expect from a photodissociation region. The most 
extreme object is MS0735, owi ng to its faint (and unc erta in) IR continuum. The larg e H2 luminosity from ofF-nuclear 
regions in the BCG NGC1275 (| Johnstone et al.ll2007[ ) led iFerland et al.l (|2008ll2009^ to propose that much of the H2 
luminosity in BCGs located in X-ray cool-core clusters can be generated by cosmic ray heating, or by non-radiative 
processes such as plasma waves. 

Brightes t cluster ga l axies are not the only galaxies to exhibit unusually large luminosities of rotational molecular 
hydrogen. lOgle et aJ] (|2007l ) find the FR II radio galaxies have strong II2 lines, but these galaxies are dissimilar to 
the BCGs in our sample. For example, 3C 326 exhibits high ionization [Ne V] and [O IV] emission, indicating AGN 
or LINER-like lines, very tiny star formation rates (< O.IM0 yr^*), and the H2 line transitions are primarily S(0) 
and S(l), indicative of cooler molecular gas than in our sample. These t ransitions are also seen in IRS mapping of 
the nearby group of galaxies, Stephan's Quintet, which exhibits bright H2 (jCluver "eFaniMoh. Sim ilarly. H2 S(0) and 
S(l) emission lines h ave been reported from IRS mapping of edge-on spiral galaxies jLaine et aLll2010t ). An archival 
study of ULIRGs by iZak amskal ()2010[ ) suggests that their H2 emission is not associated with star formation. While 
many of these studies speculate that shocks might be a source of energy (e.g. iQgle et al.ii2010i) . and might be quite 
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Fig. 11. — [Ne III] / [Ne II] ratios, indicators of the hardness of the incident UV spectrum for the ionized gas, are plotted against the 24 
lira luminosity, a surrogate for star formation rate. The ionization level of the gas shows no correlation with the mid-IR continuum. The 
neon line ratios exhibit intrinsic scatter, possibly evidence that either the ages of the youngest stars vary from system to system (which 
would mean a constant star formation rate is not a good assumption), or the mix of heating mechanisms (between star formation and 
energetic particles) differs from galaxy to galaxy. 

common, the unifying thread to all of their discussions is that the molecular hydrogen rotational lines are surprisingly 
bright and their source of energy is still unidentified. The situation is not much different here, except the BCGs tend 
to also exhibit rotational lines characteristic of warmer molecular gas than the groups or radio galaxies (S(2), S(3), to 
S(7)). 

The mid-IR luminosity is not significantly correlated with the summed luminosity of the molecular hydrogen lines, 
here represented by the sum of S(2) and S(3) lines, which were reliably detected in all 9 systems (Figure [T^). While 
we showed in § 15.21 that [Ne II] emission is correlated with dust continuum emission, here we see that dust continuum 
is not significantly correlated with molecular hydrogen emission (for fluxes, r — 0.5; for luminosities r = 0.68). We 
plot the ratio of H2 sum to mid-IR luminosity (Figure [T^b) . These ratios decrease for the systems with the highest 
mid-IR luminosities. This trend appears because the H2 luminosities are limited in range (factor of 20) while the IR 
luminosity spans a large range (> 1000). We interpret the trend to mean that the H2 heat source is more important 
and in fact dominates the IR emission features in systems with low mid-IR luminosities. 

The rather insignificant correlation between H2 and L24 utterly vanishes once JVIS0735 is omitted from the sample. 
A similar effect happens when [Ne II] fiux is compared with H2 fiux: dropping JVIS0735 from the sample causes a 
very weak (less than 2a) correlation to completely vanish. On the other hand, the luminosity of [Ne II] is correlated 
with the luminosity of the H2 S(2) -I- S(3) lines (r = 0.92). The presence or absence of 1VIS0735 has little effect on 
the inferred strong correlations between Ne II, PAHs, and mid-IR continuum flux and luminosity correlations. This 
correlation analysis suggests that while there is some relationship between the heat sources for the ionized gas and the 
dust, there appears to be a much weaker relationship between the heating sources for the molecular hydrogen and the 
dust. The ratio of molecular hydrogen to IR luminosity (Figure [T2b ) decreases with increasing IR luminosity, however, 
and suggests that the II2 heating process becomes less important to the total luminosity budget as star formation 
increases. 

The strong luminosity correlation between [Ne II] and H2 required further investigation, since the lack of correlations 
in the flux quantities suggested the luminosity correlation may be a simple "bigger is bigger" luminosity-luminosity 
comparison. (See Kennicutt (1990) for an infamous description of this type of error, involving a cigar.) Intriguingly, we 
find that the ratio of [Ne II] to L24 correlates even more strongly with the ratio of II2 S(2)-|-S(3) summed luminosity 
to L24 (r = 0.98). The best power-law fit to this relationship is H2/IR - ([Ne II]/IR)i-49±o.i2 (Figure [HI). The 
correlation of these ratios suggests that whatever process heats the molecular hydrogen is likely to be the culprit that 
boosts the forbidden line luminosity (heating the ionized gas) as well. 

To explore this idea further, we investigated how much more luminous the neon and molecular hydrogen lines are 
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Fig. 12. — Rotational H2 and infrared continuum properties compared. On left, the summed luminosity of the two most prominent 
rotational hydrogen lines (S(2), S(3)) is very weakly correlated (< 2cr significance) with the L24 as might be expected if the sample is 
(approximately) flux-limited; if plotted as a flux-flux diagram all correlation disappears. A dotted line of slope unity is shown. On right, 
the ratio of molecular hydrogen to L24 decreases with increasing L24,, which suggests that the source of energy powering H2 is unrelated 
to star formation in most of this sample. At the highest IR luminosities (A1068 and A1835), the H2 might be generated by star formation 
processes. 
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Fig. 13. — [NeII]/IR ratio correlates with H2/IR, which might be expected if the heating source that powers the molecular hydrogen well 
over that expected from star formation also elevates the [Ne II] emissivity of the ionized gas over that expected from star formation. This 
additional heating source might simultaneously explain the H2 luminosity and the excess [Ne II] luminosity. The dashed line indicates a 
slope of unity. 

than one would expect from a star forming galaxy, if the dust luminosity were a reliable indicator of the level of 
star formation. In Figure I14| we estimated the star formation rates inferred from the IR continuum (the best fits 
to the Gro ves et al. 2008 m odels, which are consistent wit h Kennicutt (1998) estimates), from the [Ne II]-|-[Ne III] 
luminosity ()Ho &: Ketoll2007D . and from the H2 luminosity (jTrever et al.ll20ldl) . For the latter estimate, we assumed 
that the Ne+ /Ne = 0.75 and Ne++/Ne = 0.15. The H2-based star formation rate in Treyer et al. (2010) relies on the 
sum of the S(0), S(l), and S(2) transitions, which were not all detected in our systems. The sums plotted are based 
only on luminosities of the detected lines. These plots show that for most of the BCGs, while [Ne II] is moderately 
over-luminous for the inferred IR-based star formation rate (a factor of 2-5 above the upper end of the scatter exhibited 
by the galaxy sample of Ho & Keto 2007, and a factor of ~ 10 over the mean), the H2 luminosity is a factor of 5-15 
over-luminous based on the IR-based star formation rates. The BCGs in A1068, A1835, and Hydra A have ratios 
typical of starbursts. 

It is interesting that the best fit for the points in Figure [13] is nearly linear. The slightly steeper than linear fit might 
be explained in the context of heating by suprathermal particles if the luminosity of the ionized gas ( [Ne II] ) is limited 
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Fig. 14. — Ratios of inferred star formation rates demonstrate that both the forbidden neon lines and the molecular rotational lines are 
being emitted in excess of what would be expected of a star-forming galaxy where the star formation rate is tracked by the IR luminosity. 
The left plot shows the ratio of the star formation rates from the sum of [Ne II] + [Ne III] (Ho & Keto 2007) compared to our estimated IR 
star formation rates, based on fits to the infrared continuum. The lower line shows the average of the Ho & Keto rate compared to the 
IR rate based on Kennicutt (1998) quantities for their sample. The higher line shows the upper limit of their scatter. Therefore [Ne II] 
in BCGs is moderately over-luminous, by a factor of ~ 3 over the mean for star-forming galaxies. The rotational transitions in H2 are 
overluminous by a factor of 5-15. In both these relations, MS0735 is an extreme example. 

by the finite column density of ionized gas, while the luminosity of the rotational line emission from the molecular gas 
is limited by the penetration depth of the suprathermal particles into the molecular gas, not the total column density 
of molecular hydrogen. 

We defer a full discussion of the excitation analysis of the individual molecular hydrogen lines to a future paper. The 
more approximate dual temperature fit that we have done here, however, shows very similar trends to those seen in 
NGC1275 filaments: the H2 rotational line intensiti es cannot be fi t by a single temperature. This trend is consistent 
with any model with a non-radiative energy source (jFerland et al.| [2008). 

In summary, as we examined correlations of continuum, PAH, and emission lines of [Ne II] and H2, and compared 
them to correlations and infrared line ratios in other types of galaxies, it emerged that a single heating process cannot 
explain the range of infrared properties we see in these BCGs. Star formation seems to play a role, albeit with 
varying levels of dominance, in producing the emission from these systems, but other processes unrelated to star 
formation must also contribute, particularly in systems with apparently low rates of star formation but high fluxes 
from rotationally-excited transitions of molecular hydrogen. 



5.4. PAHs and Dust Grain Survival and Processing 

If the dust in these BCGs spent much time in contact with the hot. X-ray emitting gas (or more generally, suprather- 
mal electrons), one might expect the dust properties, such as its size distribution or ionization fraction, to be different 
from dust that has not undergone such a traumatic experience. PAH survival alone is problematic if suprathermal 
particles alone provide heat: radiation and collisions make PAH lifetimes in the har sh environm ent of the center of a 
cool-core cluster of galaxies quite short. Using order of magnitude cross sections from' Voitl ()1992a[ ). and 0.5 keV photon 
fluxes of about 10® cm^^ s^^ we estimat e lifetimes of order one million years. The damage from particle collisions may 
be even more dire. From the analysis of IMicelotta et ahl (|2010l ). the lifetime of PAH molecules embedded in '--^ 1 keV 
gas with a density of ^ 0.1 cm""^ is limited to hundreds of years by collisions with the hot electrons and ions. Any 
processing along these lines causes the PAHs and small grains to evaporate preferentially compared to large grains. 
The presence of PAHs requires the dusty gas to be shielded from the hot gas and its radiation. 

As in other galaxies observed with the Spitzer IRS (Smith et al. 2007b; Kan eda e t al. 2008), the fluxes and lumi- 
nosities of the PAH complex at 7.7 /im and at 11.3 ^m are strongly and linearly correlated. The 7.7 and 11.3 /xm 
PAH complex luminosities are strongly correlated (r = 0.98, fluxes at r = 0.94) for all 7 systems in which both are 
detected. 2A0335 and MS0735 lack PAH 7.7 fim detections. ' 

The ratio of 7.7 to 11.3 fim PAH complexes is relatively insensitive to the sizes of the PAHs fS chutte et aHflggsl) . 
but sensitive to the ratio of ionized to neutral PAHs (Allamandola et al. 1989; Draine & Li 2007). For example, the 
ratio of the PAH complex at 7.7 to the PAH com plex at 11.3 um is lower in galaxy centers that are AGN-dominated 
compared with those which are Hll-dominated fjSmith et al.l l2007bl ) . The mean range for seven ratios of the PAH 
complexes at 7.7 and 11.3 ^m in ou r sample is 2 . 7 ib 0.2 (Figure [T5|) . The mean sample ratio is intermediate betwe en 
that of Hll/starburst galaxies (4.2, iSmith et all ()2007bf )l and diffuse Galactic emission (2-3.3, iSakon et all (|2004D 1. 
The BCGs with classic indicators of starburst activity (Abell 1068 and Abell 1835) have ratios typical of starbursts; 
the others are closer to that of Galactic ISM. In elliptical galaxies, the PAH 7.7/11.3 ratio is unusually weak (~ 1 — 2) 
(jKaneda et al, 2008) ; the lowest ratios seen in dusty ellipticals are lower than the ratios detected in our BCG sample. 
For BCGs with a low PAH 7.7 /im to 11.3 /xm ratio, the incident spectrum on the PAHs may be dominated by evolved 
stars, compared to a harder, ionizing spectrum with contributions from young massive stars. Our result suggests 
that the radiation fields in BCGs are intermediate in their hardness between those of dusty elliptical galaxies and 
star-forming galaxies, which is consistent with what one might expect in the ISM of galaxies with enormous old stellar 
populations, together with small numbers of recently formed stars. 
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Fig. 15. — The PAH complex at 7.7 fim is strongly correlated with the PAH complex at 11.3 fj,m. The mean ratio for starbursts, indicated 
with a single horizontal line at the top of the plot, is ~ 4. Diffuse Galactic ISM generates ratios from 2 — 3.3, shaded pink (light grey), and 
the ratio in dusty ellipticals ranges from 1 — 2, shaded blue (darker grey). The weighted mean for seven BCGs with PAH 7.7 fim detections 
in our sample is 2.7 ± 0.2, indicating radiation fields intermediate in hardness between dusty elliptical and normal star-forming galaxies. 



The ratio of PAH complexes at 17 and 11.3 /j,m is thought to b e regulated by PAH sizes with large PAHs contributing 
more to the longer wavelength complex (e.g.. iDraine fc Lill2007[ ). The broad PAH feature at 17 fj,m was detected with 
confidence in 4 BCGs (2A0335, A1068, Hydra A, and PKS0745), while the others have 3cr upper limits (Figure [T6l) . 
With the single exception of Abell 1068, the detections and upper limits ar e consistent with the ratio of PAH 17 //m 
to PAH 11.3 fim fluxes typical of normal star-forming galaxies 0.5) (Sm ith et al.ll2007H ). The consistency of the 
observed ratios and limits in nearly all of these BCGs compared with those with normal galaxies suggests that the 
PAH size distribution may be normal. The exception of Abell 1068, with a ratio > 2, indicates it may have deficit of 
small PAHs compared to large PAHs. This BCG also has a smaller PAH to Ltir ratio. In the BCG of Abell 1068, 
PAHs, and particularly the small PAHs, may have been preferentially destroyed by collisions with particles or photons. 

Since both [Ne III]/[Ne H] and the ratio of PAH 7.7 /zm to 11.3 /im flux are related to the hardness of radiation, 
we compared these quantities to see whether they are correlated (Figure [T7|) and whether the BCG points cover 
similar parameter space as other galaxies. We see no correlation between the PAH ratios and the neon forbidden 
line ratios, but that is not surprising since the gas producing the forbidden lines is different from the gas hosting 
the PAHs; similarly, the photons ionizing neon are not the same photons setting the ionization level of PAHs. As 
shown in Figure flTl the seven BCGs with detections in all 4 quantities exhibit ratios rather similar to SINGS galaxies 
([Smith et ;il] l2007bf). There are spiral galaxies in the SINGS sample in Figure llTl with very high [Ne III]/[Ne II] 
ratios (> 10); the explanation may be that these ga laxies (about 15% of the SINGS saniple) are dominated by ver y 
recent star formation, and therefore hotter O stars (iThoriilev et alJ [20001 : iRigbv fc Riekell200^ iSniiders et al.H 2007). 
compared to the BCGs. LINERs in the Smith et all (|2007bD sample do not exhibit [Ne III]/[Ne II] < 0.7-0.8 together 
with low PAH 7.7/11.3 /zm fractions (< 2), so a few of the BCGs in our sample with low [Ne III]/[Ne II] ratios (A1795, 
A478, and PKS0745) a lso have lower P AH 7.7/11.3 /im ratios tha n seen in t he SINGS sample, more similar to those 
of dusty ellipticals (Ka neda et al.l[2008l ) or diffuse Galactic ISM (Sa kon et al.l 2004). The two systems with the highest 
IR and PAH luminosities of the sample (Abell 1068 and Abell 1835) also exhibit the largest PAH 7.7/11.3 /im ratios, 
indicating the PAHs in the most luminous BCG systems are more ionized than those in the other BCGs. 

We arrive at the following conclusions about the PAHs and dust in these systems. 

• The presence of PAHs, and the similarity of the PAH emission ratios to those in star-forming galaxies, mean 
that these tiny grains must be protected from the ICM and shocks. PAHs are easily destroyed by ionizing UV 
and X-rays, by collisions with hot thermal particles, and by shocks. 
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Fig. 16. — The luminosity ratios and upper limits for the PAH complexes at 17 and 11.3 fim. The line indicates the mean value of the 
ratio between the 17 and 11.3 fim PAH complexes (0.5) for normal star-forming galaxies (Smith et al. 2007^). Only A1068 has a ratio that 
is distinctly atypical, suggesting some processing of the smaller PAHs have occurred in that system. 

• The emissivities of the PAHs and the ionized gas are correlated with the mid-IR continuum emitted by dust 
grains. Therefore there is a common source of heat for these components, consistent with being star formation. 
However, the excess luminosity of [Ne H] and H2 in the less luminous systems suggests that another component 
is contributing to the heating of the ionized gas in addition to star formation. 



• The emissivity of the PAHs is not correlated with the stellar (6 /im) IR continuum, 
the stars were the main agent heating the dust. 



Correlation is expected if 



We speculate that the evolved stars in the BCGs are the main production source of dust. The dust, except in a 
couple of cases, shows little indication that it may have been processed by hot, X-ray emitting plasma or AGN. The 
source of the cool BCG gas may be the hot intracluster medium, but it seems unlikely that the dust came with the 
gas. Dust, however, is an extremely effective coolant, and if mixed with the hot gas, could precipitate the cooling 
necessary to fuel the formation of dusty molecular star-forming clouds. 



5.5. AGN Contributions 

Synchrotron radiation from a jet could contribute some of the infrared continuum. iClearv et al.l ()2007| ) find in their 
survey of an unbiased set of SCR and quasars that only the quasars have a nonthermal contribution at 15 /xm of 
> 20%. Abo ut half of the 15 -Mm emission from nucleus of M87 is synchrotron emission, based on analysis of its IRS 
spectrum (Bu son et al.l 120091 ) . Of the 9 sources in our sample, only Hydra A has a radio synchrotron source that is 
luminous enough to contribute more than ^ 10% to the IRS continuum at 24 /im, if its power law extends unbroken 
into the infrared. 

High ionization lines in the infrared could indicate a buried AGN or the presence of very hot star s. [Ne V] 14 .3 iim 
and 24.3 fiui, with an ionization potential of 97.1 eV, would be an unambiguous AGN indicator (|Voitill992bD . We 
did not detect this emission line in any of our sources at a level > 10""'^'' W m~^), although a a weak [Ne V]14.3 /xm 
feature is detected in Abell 1835 and Hydra A, and a suspicious (possibly spurious) feat ure appear s near 24 /im in the 
spectrum of Abell 1068, which is known to have some AGN contribution in the infrared (jQuillen et al. 2008). The lack 
of [Ne V] does not rule out some AGN contribution, since this line can be faint compared to [Ne II]. [O IV] at 25.9 fj,m 
would also be an unambiguous sign of AGN photoionization, since O HI has an ionization potential just above that 
of He I (to He II), 54.4 eV, and is therefore rare in regions ionized by stars. However, [O IV] 25.9 fj,m is blended with 
[Fe II] at 26 /im. We detected this blend unambiguously only in PKS0745-19 and Hydra A. Puzzlingly (if this blend 
is an indication of the presence of O IV), the [Ne III]/[Ne II] ratio in PKS0745-19 is one of the lowest in the sample. 
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their Table 4) with black points, elliptical galaxies (Kaneda et al. 2008) with blue filled circles inside triangles, and our sample (red, filled 
square points). The [Ne III]/[Ne II] line ratios of BCGs arc similar to those seen in star-forming galaxies with low [Ne III]/[Ne II] ratios, 
LINERS, and dusty cUipticals. Three of the red (filled square) points (A478, A1795, PKS0745) have somewhat lower PAH 7.7 to 11.3 
^tm ratios than seen in SINGS galaxies with similarly low [Ne III]/[Ne II] ratios, but are rather similar to the ratios seen in ellipticals in 
IKaneda etHI l l200l) . 2A0335 and MS0735 are not plotted here. 

The [O IV]- [Fe II] blend docs not appear at all in the spectrum of Abell 1068, however, nor does [Ne V] 14.3 /im. The 
lack of these features casts further doubt on the candidate [Ne V] feature in that spectrum. 

The presence of high ionization lines such as [S IV] (10.51 fim) and [Ne III] indicate the presence of young hot stars 
in these galaxies. The ionization parameters consistent with the [Nelll] / [Nell] ratios in the ionized gas are similar to 
those seen in star-forming galaxies. The energy required to ionize neon once to Ne II is 21.6 eV, compared to 41 eV 
for reaching Ne III. [Ne II] was easily detected in every galaxy in our sample, and only ]V[S0735-t-74 lacked detectable 
[Ne III]. Although S III to S IV has an ionization potential (34.8 eV) similar to Ne II to Ne III, it has a very low 
photoionization cross-section. So a significant detection of [S IV] not only indicates hot stars, but a high density of 
them. In contrast to [Ne III], [S IV] was not detected in any of the systems. PAHFIT resuhs for 2A0335-I-096, Abell 
1835, Abell 1068, and Hydra A indicate faint formal detections at 3 — 7a, but inspection of the fits and data causes 
us to regard these detections as extremely marginal. 

Based on the lack of [Ne V] emission lines from the gas, we have no conclusive evidence in favor of the gas being 
photoionized by an AGIV, consistent with concl usions based on spatially-resolved optical emission-line studies of these 
and similar galaxies (e.g.. iHeckman et"allll989| ). 

The detection of spatially extended UV continuum in a number of case s is unambiguous evidence of the importance 
of recent star forrn ation ove r AGN contrib utions i n this sample of BCGs (|]V[cNama,ra fc 0'Connelliri993t iMartel et alJ 
[200llO'Dea et al.ir2 004: Hic ks fc IVIushotz kv 2005: Hicks eTal|[20lol: iDonahue et al.ll2010[ ). Sincelack of evidence is not 
the same as evidence of lack, we keep in mind that some of the spectra may have contributions from a low-luminosity 
AGN since none of these spectra exclude the nucleus. However, Occam's Razor prefers the simplest explanation for 
heating these systems, and as such, no need for AGN excitation or heating is required by these observations. 

6. CONCLUSIONS 

In our Spitzer IRS spectroscopic study of 9 brightest cluster galaxies residing in cool-core X-ray clusters, we have 
detected very bright molecular hydrogen rotational transitions, PAH features, and forbidden lines from ionized gas 
(Ne II, Ne HI), in addition to dust continuum at 15-25 / zm. These galaxies we re known previously to have prominent 
optical forbidden and Balmer emission line nebulae (e.g. IHeckman et al.lll989l) . Photometric IVIIPS and IRAC studies 
of similar BCG s have shown prom inent mid-infrared dust emission, together with UV and blue signatures of star 
formation (e.g. lO'Dea et al.l 120081) . The low resolution Spitzer 5-25 fj,m IRS spectra reveal that BCGs host PAHs, 
and the ratios of these features indicate that the PAHs in such BCGs are similar to PAHs in other types of galaxies. 
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The emission from PAHs and dust are highly correlated, nearly linearly, as might be expected if the PAHs and dust 
are heated by star formation. The ratios of the 11.3 iim PAH luminosities to total IR luminosities are similar to or 
slightly lower than those of normal star- forming galaxies, with the exception of MS0735 whose PAH 11.3 fj,m feature 
is well-detected, and its IR continuum is very weak. Fits of simulated starburst models published in Groves et al. 
(2008) show that the star formation rates inferred from the 5-25 /im spectra are consistent with star formation rates 
inferred from our 70 /im MIPS photometry using relations from Kennicutt (1998) and Calzetti et al (2010). These 
simulated starburst SEDs also provide reasonable fits to the PAH features in the spectra, again consistent with the 
strong, nearly linear correlation between PAH luminosity and dust continuum luminosities. 

The luminosity of the warm ionized gas, traced by the [Ne II] 12.6 fim emission line, is also correlated with the dust 
luminosity, so it is likely that the star formation powering the dust luminosity is also contributing to the heating of 
the warm ISM. However, the relation is distinctly non-linear, in the sense that the systems with lower IR luminosities 
have larger [Ne II] /IR ratios. All of the systems are over- luminous in [Ne II] compared with PAH or dust emission, as 
most l ie at or above the top end of the scatter of this luminosity relation as seen in star-forming galaxies ()Smith et al.l 
I2007bfl . 

Even more strikingly, the molecular hydrogen luminosities of BCGs are very high compared to that expected from 
star-forming galaxies of similar infrared luminosities. The H2 luminosities are only weakly correlated with the mid IR 
or PAH luminosities, suggesting that the H2 emission line power source is nearly independent of that powering the 
dust. The strong correlation of the molecular hydrogen line luminosity scaled by L24 with [Ne II] scaled in the same 
way suggests that the warm gas has a second heat source that may be related to the primary power source for the 
molecular hydrogen emission l ines, a scenar io consiste nt with excess heatin g by energetic particles and / or conduction 
from the hot intracluster gas (jSparks et al.| [T989: Fcrl and et aLll2008l[2009l ). 

The existence of dusty gas and, in particular, PAHs, along with the fact that PAH feature ratios are similar to those 
found in star-forming galaxies, suggests that (a) the dust and PAHs are similar in ionization and size distributions as 
in spiral galaxies and that (b) the PAHs and dust are shielded from the destructive X-ray radiation and fast-moving 
thermal particles in the ICM. On the other hand, the ionized and molecular gas may indeed receive a noticeable 
energy dose from thermal particles from the hot ICM, boosting the forbidden line radiation from the ionized gas and 
the rotational line emission from mo lecular hydrogen. 

The fits with starburst SEDs from ' Groves et al.r(l2008| ) give a consistent story to our empirical comparisons: we can 
simultaneously fit the mid-IR continuum and PAH features and some of the ISM emission lines, and the fits give SFRs 
which are consistent with Kennicutt (1998) relations for the 24-70 fim photometry measured with Spitzer MIPS and 
with Calzetti et al. (2010) relations for 70 fj,m. However, the models do not reliably predict the neon line spectrum, 
which we suggest may be partially generated by the same physical mechanism that generates the majority of the 
molecular H2 emission. 

Furthermore, our results are in agreement with the iJohnstone et al.l (|2007t ) analysis of IRS spectra of NGC1275 
and NGC4696. Our sample further demonstrates that there is no correlation between the star-formation signatures 
in these galaxies and the strength of the molecular hydrogen luminosities, and that molecular hydrogen luminosities 
are more extreme than produced by typical star formation-related processes. The neon forbidden lines and rotational 
lines of molecular hydrogen can be enhanced without significantly modifying the luminosities of dust if the ionized gas 
and some of the molecular gas is heated by non-radiative agents that can penetrate the molecular gas clouds, such as 
cosmic rays, hot electrons, MHD waves. 
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